Random matrix theory for underwater sound propagation 
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Ocean acoustic propagation can be formulated as a wave guide with a weakly random medium 
generating multiple scattering. Twenty years ago, this was recognized as a quantum chaos problem, 
and yet random matrix theory, one pillar of quantum or wave chaos studies, has never been intro- 
duced into the subject. The modes of the wave guide provide a representation for the propagation, 
which in the parabolic approximation is unitary. Scattering induced by the ocean's internal waves 
leads to a power-law random banded unitary matrix ensemble for long-range deep ocean acoustic 
propagation. The ensemble has similarities, but differs, from those introduced for studying the An- 
derson metal-insulator transition. The resulting long-range propagation ensemble statistics agree 
well with those of full wave propagation using the parabolic equation. 

PACS numbers: 43.20Bp, 05.45.Mt, 43.30.Ft, 43.20.Bi 



Underwater sound provides a means of remote sens- 
ing the ocean interior and monitoring global climate 
change [1, 2] amongst other motivations. Beginning in 
the latter part of the 1980's decade it was realized that 
the ray dynamics underlying ocean acoustic propagation 
in many contexts was chaotic [3-6] and therefore the field 
serves as a domain with several unique features for stud- 
ies of wave/quantum chaos [7-9]. To date, one of the the- 
oretical foundations of quantum chaos, random matrix 
theory (RMT) [10-12], has not been applied previously 
in the ocean acoustic propagation context. This is in 
stark contrast to other linear acoustics fields, where ran- 
dom matrix theory has been of growing usefulness since 
its introduction over twenty years ago [13, 14]. RMT 
has the potential to improve the understanding of many 
of the observed statistical behaviors and connect them 
to parameters of the ocean environment, and give a new 
more efficient method of making simulations. Conversely, 
unique features of ocean acoustics have the potential to 
influence the future development of RMT through the 
introduction of new models. 

The earliest RMT application to elastodynamics [13] 
had a direct experimental connection to the classical, 
structureless Gaussian/circular ensembles of Wigner and 
Dyson [12]. Ocean acoustic propagation cannot be repre- 
sented in this way. The purpose of this letter therefore is 
to step through the construction of a structured random 
matrix ensemble appropriate for ocean acoustic propaga- 
tion in as straightforward a physical context as possible. 
For simplicity, long range propagation at a fixed, low 
angular frequency u is considered. The sound channel 
creates a vertical wave guide [15] in which the ocean's 
internal waves [16] generate multiple scatterings (intro- 
ducing chaos). Horizontal out-of-plane scattering can be 
neglected, as can absorption, surface, and bottom inter- 
actions. Furthermore, the dominant effect of internal 
waves is small angle, forward scattering. Incorporation 
of these approximations into the appropriate scalar wave 



equation leads to a paraxial optical (parabolic) equation. 
It has a direct analogy to the Schrodinger equation, mak- 
ing the wave and quantum chaos connection even closer. 
In this analogy, the wave vector inverse k^ 1 plays the 
role of a pseudo-Planck's constant and the propagation 
range r the role of a pseudo-time, i.e. 
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where the real part of Sfr(z,r) multiplied by a traveling 
phase and r -1 / 2 is the sound pressure amplitude. The 
"potential" is 
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Co is a reference sound speed « 1.5km/s, and k^ 1 = 
co/uj. Vo(z) typically has a minimum at approximately 
1 km depth and is modeled with a Munk profile [15]. 
It possesses an exponential form near the surface due 
to temperature decrease and a linear dependence deep 
underwater due to pressure increase. Multiple scattering 
is induced by the internal wave fluctuations in the range 
and depth dependences Vi(z,r) and the form used here 
is that of Ref. [17]. More details and further references 
can be found in the review paper [7]. Note that our 
calculations have no cutoffs in mode number for the ocean 
surface or bottom. 

Wave field propagation for a parabolic equation is uni- 
tary and the modes provide a physically relevant basis 
with which to study the effects of scattering [18, 19]. Ex- 
tensions to adiabatically denned modes could also be in- 
corporated, but here the modes are independent of range. 
The unitary propagation matrix U expressed in the mode 
basis is the natural vehicle for which to construct the 
random matrix ensemble. The problem is to identify and 
incorporate into the ensemble's structure all information 
that survives long-range propagation (up to several thou- 
sand kilometers) and no more. 
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Without scattering from internal waves, or indeed any 
plausible scattering mechanism, the unitary propagation 
matrix would be a diagonal matrix, U d%ag = A, and ac- 
cumulate a phase proportional to the range propagated, 

AmnO) = 5 mn exp(-ik E m r) , (3) 

where E m is the energy of the m th mode of the vertical 
waveguide. To propagate an initial wave field just re- 
quires its initial modal decomposition and these phases. 

Perturbatively to first order, internal waves generate a 
propagator 

UxiA^I-iekoJ dr'V^J , (4) 

which requires restoration of unitarity (Vj is the operator 
corresponding to Vi(z, r) in the interaction picture). One 
technique applies the Cayley transform of a Hermitian 
operator A 

U = A(I + ieA)-\l - ieA) , A = / dr'Vj , (5) 
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which is quite useful because the internal waves generate 
a weak scattering locally and evaluating the operator in- 
verse is straightforward. As a consequence, it is natural 
to construct a fixed-range building block unitary prop- 
agator Ub = U(r = 77, ) for ranges long enough for the 
internal waves to generate sufficiently random behaviors, 
but short enough that perturbation theory is still a viable 
approximation. Transfer matrix building blocks are used 
in Ref. [20] . Long range propagation follows by multiply- 
ing the requisite number of independently drawn mem- 
bers of a Ub ensemble to arrive at an ensemble of U for 
the full range [i.e. U(r = £rb) = Yl^ Ubj]- 

Hence, ideally there are three statistical properties of 
interest to determine a practical value for the range rv 
i) correlations between a matrix element of Ub for a given 
block and its adjacent block should be small, ii) the 
phases of each Ub^mn should be largely randomized, and 
iii) ideally, dynamical correlations between neighboring 
matrix elements should be minimal. For extremely short 
range propagation, there is little scattering, U is nearly 
diagonal, and the matrix elements are highly correlated 
with little randomness. As the propagation range in- 
creases, there comes a point at which the phases of the Ub 
matrix elements become more or less randomized. Fig- 
ure 1 shows the randomization is well underway by 50 
km at 75 Hz (used throughout). At this range the cor- 
relations between Ub,mn for adjacent blocks have fallen 
to roughly 10% or less. Conveniently, this is roughly the 
range at which wave energy has cycled once from the 
upper to lower turning point and back again. As the 
strongest perturbations are close to the surface, choosing 
this range is consistent with accounting for one strong 
perturbation cycle. Other than effects due to the unitar- 
ity constraints, the matrix elements are beginning to be- 
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FIG. 1. The phases of the matrix elements of A~ 1/2 UA~ 1/2 
for propagation through a single internal wave field. The top 
panel is for 1 km, and the bottom panel 50 km. The nonran- 
dom behaviors disappear with increasing range. 



have like zero-centered, complex, Gaussian random vari- 
ables. Thus, the A mn are taken as independent, com- 
plex, Gaussian random variables; the diagonal elements 
are real. The dynamical correlations amongst neighbor- 
ing matrix elements are not so small, especially along 
diagonals of U^ 'dynamical' means correlations in addi- 
tion to those induced by unitarity. Nevertheless, as a 
starting point, these correlations are ignored. If later, 
it is found that the ensemble is in some way deficient, 
one could revisit their incoporation. That leaves just the 
variance determination for each matrix element. 

One might anticipate that the mode-mixing from scat- 
tering cannot involve modes separated too greatly in 
mode number rendering the matrix banded, at least ini- 
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FIG. 2. [/for propagation through a single internal wave field 
is illustrated as a color plot of the magnitudes of the matrix 
elements |C/ mn |. The top panel is for 50 km propagation, and 
the bottom panel for 1000 km. 



tially in the propagation. Figure 2 illustrates the band- 
edness of U due to the internal wave field. Thus, the vari- 
ance of the A mn decreases rapidly away from the diag- 
onal. The matrix possesses a rough translational invari- 
ance along the diagonal and the matrix element variance 
is mostly a function of \n — m\. Figure 3 displays the 
behavior of the band, which has a power law decrease of 
exponent ^1.3 for the standard deviation up to a shoul- 
der near \n — m\ ^50, where it suddenly drops off much 
more sharply. Further investigation of the shoulder is 
needed to ensure that it is not an artifact of the internal 
waves construction method. 

The ensemble of Eq. (5) with the variances of Fig. 3 
has three features distinguishing it from the power-law 
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FIG. 3. The variance of the off-diagonal matrix elements of 
U as a function of distance from the matrix diagonal. The 
solid line is the result of a perturbation theory not explicitly 
discussed in the text. The dashed line indicates the power- law 
decay exponent. 

random banded matrices introduced for investigations of 
the Anderson metal-insulator transition [21]. First, the 
ensemble is unitary; see Ref. [22] for introduction and 
discussion of a unitary ensemble. Second the diagonal 
elements retain a deterministic function and importance 
depending on the value of e. Finally, there is a shoulder 
in the band width. The band width exponent is roughly 
independent of wave vector and just a bit greater than 
unity, which makes this ensemble most similar to those 
possessing localization and super-diffusive wave packet 
spreading at short times [21]. Thus, the ensemble is con- 
sistent with localization in mode number, which allows 
for the well known possibility of isolating and measuring 
early arrival structures [23, 24]. Were the internal waves 
or other scattering mechanisms in the ocean of a differ- 
ent character with a slower decay away from the diagonal, 
that would not have been the case. The bandwidth de- 
cay exponent is little changed by multiplication of [/&, at 
sufficiently long propagation range the width will grow 
sufficiently to make wave energy hit the ocean's surface 
and bottom, which will then strip out energy. 

It is a somewhat brash assumption that dynamical cor- 
relations of neighboring matrix elements Uh.nm can be 
completely ignored. Nevertheless, it is seen in Fig. 4 
that a typical unitary banded random ensemble matrix U 
obtained as a product over building blocks and a U con- 
structed by propagation through a single realization of an 
internal wave field using the parabolic equation give the 
same general appearance, including band width, when 
propagated out to 1000 km. From this figure at least, 
there is no indication that this correlation information 
survives intact in propagation to long ranges. Further 
details are to be published [25]. 

In summary, for the first time a random matrix en- 
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FIG. 4. Illustration of U as a color plot of the logarithm of 
|C/ mn |. The upper panel results from a single ensemble real- 
ization, and the lower panel results from parabolic equation 
propagation through a single internal wave field. The propa- 
gation range is r = 1000 km. 

semble has been introduced for ocean acoustic propaga- 
tion. The ensemble is constructed for a building block 
unitary propagator, and further propagation is by ma- 
trix multiplication of independently chosen The ma- 
trix is banded in mode number with an approximate 
power law whose exponent firmly places the ensemble 
in the class of those that exhibit localization and super- 
diffusion. Unlike the ensembles introduced to study the 
Anderson metal-insulator transition [21], the ensemble 
here is unitary, has deterministic structure for the diag- 
onal elements, and there is a cut off /shoulder in the band 
width where the matrix elements fall off precipitously. 
The relatively slow decay of a power law cutoff has ram- 
ifications for convergence of calculations in determining 



where one can make numerical approximations such as 
matrix truncations. The ensemble construction is much 
faster than full wave propagation and for statistical pur- 
poses can be used to simulate the dynamics without hav- 
ing to explicitly construct the potential due to internal 
waves. There have been a number of investigations of 
power-law banded random matrices in multiple contexts, 
and there exists some analytic analysis [26] that may in- 
form ocean acoustics applications. It would be interest- 
ing to consider shallower water and derive ensembles with 
absorption and surface/bottom scattering. 
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